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ABSTRACT 

The central image of a strongly lensed background source places constraints on the foreground 
lens galaxy’s inner mass profile slope, core radius and mass of its nuclear supermassive black hole. 
Using high-resolution long-baseline Atacama Large Millimeter/submillimeter Array (ALMA) observa¬ 
tions and archival Hubble Space Telescope ( HST ) imaging, we model the gravitational lens H-ATLAS 
J090311.6+003906 (also known as SDP.81) and search for the demagnified central image. There is 
central continuum emission from the lens galaxy’s active galactic nucleus (AGN) but no evidence of 
the central lensed image in any molecular line. We use the CO maps to determine the flux limit of the 
central image excluding the AGN continuum. We predict the flux density of the central image and use 
the limits from the ALMA data to constrain the innermost mass distribution of the lens. For a power- 
law profile with a core radius of 0.15" measured from HST photometry of the lens galaxy assuming 
that the central flux is attributed to the AGN, we find that a black hole mass of F^Mbh/Mq) > 8.5 
is preferred. Deeper observations with a detection of the central image will significantly improve the 
constraints of the innermost mass distribution of the lens galaxy. 

Subject headings: gravitational lensing: strong 

1. INTRODUCTION 

Galaxies with bulges harbor supermassive black holes 
(SMBH) at their centers. The SMBH mass is correlated 
with physical properties of the bulge (e.g, luminosity, 
stellar velocity dispersion) despite the bulge extending 
beyo nd the black hole’s dynamical sphere of inf l uence 
(e.g.. IMa eorrian et all I199 8|: |Ferrarese fe Merritt! 120001 : 

IGebhardt et al.l 120001 : iGultekin et all 120091) . These sur¬ 

prising correlations sugg est a coev o lution of the SMBH 
and its host galaxy le.g.. lKormendv fe 1 ldl201 itl . Deter¬ 
mining the origin of such correlations is important for 
understanding how galaxies form and evolve. 

Direct measurements of black hole masses based on 
stellar dynamics, gas dynamics or maser dynamics are 
restricted to gal axies within ~ 15 0 Mpc (see, e.g., Ta¬ 
bles 2 and 3 of iKormendv fe Hoi 120131 and references 
there in). Reverberation mapping dBlandford fe McKcd 
Il982f) can measure black hole masses of active galactic 
nuclei (AGNs) at distances up to ~ 1 Gpc, especially for 
the m ost massive black holes of ~ 1O 9 M 0 ('e.g. JPetersonl 
IM1 . However, reverberation mapping is only possi¬ 
ble for galaxies with bright AGN emission, complicat¬ 
ing measurements of the velocity dispersion and stellar 
mass of the host galaxy’s bulge. Beyond 2 ~ 0.4, one 
must rely on assumed scaling relations between SMBH 
mass and some observable property of the AGN (e.g., 
luminosity), as direct measurements are not possible at 
these cosmological distances. One way to independently 
measure SMBH masses at these distances and beyond is 
through strong gravitational lensing. 

Strong gravitational lensing occurs when a massive 
foreground object is located close to the line of sight to a 
background source. By fitting a mass model to the mul- 
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tiple lensed images of the background source, properties 
of the mass distribution of the foreground lens galaxy 
can be inferred. These constraints are strongest within 
the galaxy’s Einstein radius, which is on the scales of a 
few kpc. Lensing theory predicts that for a non-singular 
mass distribution, there should be an additional highly- 
demagnified image very close to the center of the lens. 
The brightness of this image is highly sensitive to the 
central mass distribution of the lens on much smaller 
scales (~ 100 pc), with more concentrated mass distri¬ 
butions producing greater demagnification. Detection of 
these central images is extremely challenging due to their 
low brightness and the fact that they are embedded in 
emission from the lens galaxy. The best prospects for 
detection are in radio observations, where lens galax¬ 
ies generally have very weak emission. Indeed, the only 
conf irmed detectio n of a central image is of a radio lens 
(jWinn et al.l I2004T) . However, even a non-detection can 
place constraints on the inner mass distribution (e.g., 
IZhang et al.ll2007h . 

Recent studies reveal a population of lensed dusty star¬ 
forming galaxies at z ~ 2 — 3 (e.g.. INeerello et al.ll2 010t 
Bussmann et al.l 120131 : iHezavch et al.l 120131 : lYieira et al l 
2013T) . Due to the steepness of their luminosity func¬ 
tion at the bright end, these galaxies benefit greatly 
from lensing magnification. Additionally, the brightness 
of these galaxies is relatively unaffected by distance be¬ 
yond z > 1 due to the negative K-correction, where 
the rising SED from the dust emission compensates the 
dimming due to increasing cosmological distance. As 
a result, these lensed dusty star-forming galaxies are 
among t he brightest sources in w ide-area submillimeter 
surveys dHezaveh fc Holdeil 1201 Ilf . Their large submil¬ 
limeter fluxes make them ideal targets to search for cen¬ 
tral images. With the high sensitivity and resolution 
of the A tacama Large M illimeter/submillimeter Array 
(ALMA! . iHezaveh et all (1201511 predict that observations 
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of these lenses could detect their central images and con¬ 
strain the size of the lens galaxy’s core, mass profile slope, 
and mass of its SMBH. 

In this paper, we use the first long-baseline ALMA 
observations of a strong gravitational lens, H-ATLAS 
J090311.6+003906 (hereafter SDP.81), to place an up¬ 
per limit on the flux of the central image and constrain 
the properties of the innermost regions of the lens galaxy. 
SDP.81 is a massive elliptical galaxy at z = 0.2999 lens- 
ing a background dusty star -forming galaxy at z = 3.04 2 
into four multiple i mages dN egrello et all 1201(1 120141 b 
It was identified by iNegrello et all (|2010f) as one of the 
brightest sources in the Science Demonstration Phase 
(SDP) of the Herschel Astrop hvsical Terahertz Large 
Area Survey ('Ll-ATLAS: lEales et al.ll2010lb 

This paper is organized as follows. In Section [2J we 
summarize the data used in this analysis. Our lens mod¬ 
eling procedure is described in Section [3[ We present our 
results in Section 3] and summarize in Section [5] We as¬ 
sume f lm = 0.3, Da = 0.7, and Ho = 70kms -1 Mpc -1 . 
All quantities are given in h^o units. At z = 0.2999, the 
angular scale is 1" = 4.45 kpc. 

2. DATA 

2.1. ALMA Data 

ALMA Science Verification observations of SDP.81 
were taken in October 2014 as part of the 2014 ALMA 
Long Baseli ne Campaign. Deta i ls of t he observations are 
described in [Partnership et all (120151b We use archival 
images of the CO lines and the Bands 4 (151 GHz) 
and 6+7 (268 GHz) continuum, CLEANed with the 
lOOO/cAitu tapering. For the Bands 6 (236 GHz) and 
7 (290 GHz) continuum images, we created lOOOkXuv 
tapered images from the archival calibrated uv data in 
the same manner as the other lOOOfcAuu tapered images. 
The images are shown in Figure [l] 

The beam sizes and position angles (measured East 
of North) are as follows: 0.182" x 0.143", 57.3° for CO 
J=5-4, 0.199" x 0.138", 46.8° for CO J= 8-7, 0.202" x 
0.131", 42.4° for CO J=10-9, 0.158" x 0.139", 61.0° for 
the Band 4 continuum, 0.192" x 0.133", 46.5° for the uv- 
tapered Band 6 continuum, and 0.200" x 0.125", 42.5° for 
the uv-tapered Band 7 continuum. The root mean square 
(rms) noise level of each map is 0.0166 Jy beam -1 km s -1 
for CO J=5-4, 0.0194 Jy beam -1 km s -1 for CO J= 8-7, 
0.0214 Jy beam -1 km s -1 for CO J=10-9, 0.0115 mJy 
beam -1 for the Band 4 continuum, 0.0200 mJy beam -1 
for the uv-tapered Band 6 continuum, and 0.0193 mJy 
beam -1 for the uv-tapered Band 7 continuum. 

There is compact emission from the center of the 
system , which was also identified by iPartnershin et aP 
•: 2015 :. We distinguish between a central image and 
emission from the lens galaxy by compar ing its SEP to 
that of the primary le nsed images (e.g., iMcKean et al.l 
[20051 IMore eTaIll2008lb Its flat SED (see Section [3]) in¬ 
dicates that it is not the central image, but rather low- 
level AGN emission from the lens galaxy (Figure [2J . The 
central image is not detected in any of the continuum- 
subtracted line emissions, indicating that it is highly de- 
magnified. Being a quad lens with four bright images (A, 
B, C and D in Figurejl]) of the background source, this is 
expected for SDP.81 since quads generally have more de- 
magnified central images than doubles dMao et al.ll2Q01k 


lKeetonl[200l . 

2.2. HST Data 

We use archival Hubble Space Telescope ( HST) Wide 
Field Camera 3 (WFC3) imaging of SDP.81 to model the 
light distribution of the lens galaxy and constrain its core 
radius. The observations were taken in April 2011 (pro- 
p osal #12194, P I: Negr ello). These data were presented 
in Neg rello et al.l (120141 a nd were used to model the lens 
system dDve et al.l 12014 b The observations consist of 
712 seconds of exposure time in F110W and 4418 sec¬ 
onds of exposure ti me in F160W. The da ta are reduced 
using DrizzlePac dGonzaga et al .1 [201 2fPl with resam¬ 
pling to a 0.065"/pixcl scale. The data are shown in the 
top panels of Figure OH 

3. LENS MODELING 

Our lens modeling is performed with Glee, a software 
developed by S. H. Suyu and A. Halk ola (|Suvu fe Halkolal 
120101 : Isiivu et, al.ll2012h . In Section 13.11 we identify the 
positions and fluxes of the multiple images of the lensed 
background source that are used as constraints for model¬ 
ing our lens mass distribution. In Section [3^1 we further 
use the HST imaging to place constraints on the inner 
mass distribution of the galaxy. Lensing mass distribu¬ 
tions are parametrized profiles, and we explore two forms 
of mass model in Section [+3l (1) a total mass distribu¬ 
tion that follows a cored power law distribution, and (2) 
a total mass distribution composed of baryons and dark 
matter. Model parameters of the lens are constrained 
through Markov Chain Monte Carlo (MCMC) sampling. 

3.1. Constraints From the Multiple Images 

We visually identify three distinct features in the 
lensed source from a combination of the high-resolution 
Band 6 and Band 7 continuum, the lower-resolution 
Band 4 continuum, and the velocity-resolved Band 4 
data. The locations of the identified features are shown 
as open symbols in Figure [4] on top of the CLEANed 
high-resolution Band 7 image. Feature 1 is identified 
as a quadruply-imaged compact star forming clumps in 
the high-resolution Band 7 continuum image. Features 
2 and 3 are identified as doubly-imaged emission regions 
in the Band 4 continuum image. These three features 
correspond to the thre e distinct c lumps identified by 
iRvbak et al.l (|2015h and iDve et al.1 (|2015| b one of which 
is inside the tangential caustic and two of which are out¬ 
side. The image positions used as constraints on the lens 
model are given in Table [I] 

The flux ratios of individual features are difficult to 
measure since these features are close in proximity, and 
it is not straightforward to separate these features for 
flux measurements as the fluxes depends sensitively on 
the choices of regions for separating these features. We 
therefore measure fluxes of the entirety of each lensed 
image (FigureQ] upper left panel). We designate the four 
images A (left arc), B (lower left arc), C (upper left arc), 
and D (right counterimage). We designate an additional 
region E, where the lens center is located and where we 
expect to detect a central image if possible. We use the 
Band 4 continuum, uv-tapered Band 6 continuum, and 

3 http://drizzlepac.stsci.edu/ 
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Fig. 1.— ALMA images of SDP.81. Shown are the uv-tapered Band 7 continuum (top left), uv-tapered Band 6 continuum (top center), 
Band 4 continuum (top right), CO J=10-9 (bottom left), CO J= 8-7 (bottom center), and CO J=5-4 (bottom right). The beam size is the 
white ellipse in the bottom left of each panel. The red regions in the Band 7 continuum image (labeled A-E in order of decreasing flux 
density) are used to compute the flux densities of the images. All images are 4.2" on a side. 


TABLE 1 

Image Position Constraints 


Image ID 

a (J2000) a 

<5 (J2000) 

a 

Uncertainty ( 7/ ) 

At 

09:03:11.673 

(-1.505) 

+00:39:06.40 (- 

-0.130) 

0.03 

B1 

09:03:11.628 

(-0.830) 

+00:39:05.34 (- 

-1.190) 

0.03 

Cl 

09:03:11.645 

(-1.085) 

+00:39:07.54 

(1.010) 

0.06 

C2 

09:03:11.629 

(-0.845) 

+00:39:07.64 

(1.110) 

0.03 

C3 

09:03:11.619 

(-0.695) 

+00:39:07.67 

(1.140) 

0.03 

D1 

09:03:11.446 

(1.900) 

+00:39:06.74 

(0.210) 

0.03 

D2 

09:03:11.444 

(1.930) 

+00:39:06.46 (- 

-0.070) 

0.03 

D3 

09:03:11.445 

(1.915) 

+00:39:06.18 (- 

-0.350) 

0.06 


a Values in parentheses indicate offset in arcseconds relative to the mass centroid prior. 


uv-tapered Band 7 continuum images to estimate flux 
ratios. Uncertainties on the fluxes in these regions are 
the background rms multiplied by the square root of the 
number of beams in the region. The fluxes of the regions 
are given in Table [2] 

We note that images A and B only contain flux from 
the portion of the source located within the tangential 
caustic. Additionally, since the background source is 
spatially extended and the magnifications of images A, B 
and C are high, we do not use the fluxes of these images as 
constraints due to t he high u ncert ainty from differential 
magnification fe.g.. iHezaveh et al.H2012l) . which we esti¬ 
mate to be ~ 30-50% from simulations. We use only the 
positions of the three identified features as constraints on 
the lens model. We then use the flux of D and the flux 


limit of E in the CO images for constraining the inner lens 
mass distribution. The use of continuum-subtracted CO 
fluxes ensure that they originate from the lensed back¬ 
ground source, excluding other contaminating emissions 
such as those from AGN in the foreground lens galaxy. 


3.2. Light Fitting to HST Data 

In order to model the inner mass distribution of the 
lens galaxy, we first fit light profiles to an 87 x 87 pixel 
(5.655" x 5.655") region of the WFC3 data, excluding an 
annular region containing the lensed images (Figure [3]). 
A total of 3575 pixels are used for the fitting. We fit the 
lens light with the sum of two cored elliptical power-law 
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TABLE 2 
Image Fluxes 


Image 

CO J= 5-4(Jy kms- 1 ) 

CO J=8-7(Jy km s" 1 ) 

CO J=10-9(Jy km s J ) 

Band 4 (mjy) 

Band 6 (mJy) 

Band 7 (mjy) 

A 

2.61 ±0.13 

2.58 ±0.16 

0.97 ±0.14 

1.33 ±0.08 

7.54 ±0.14 

10.8 ±0.14 

B 

2.16 ±0.12 

2.38 ±0.14 

0.80 ±0.12 

1.11 ±0.08 

5.78 ±0.13 

8.20 ±0.13 

C 

1.81 ±0.11 

2.03 ±0.13 

0.55 ±0.11 

1.12 ±0.07 

5.32 ±0.12 

8.09 ±0.12 

D 

1.45 ±0.11 

1.92 ±0.13 

0.59 ±0.11 

0.99 ±0.07 

5.04 ±0.12 

7.64 ±0.11 

E 

0.03 ±0.03 

0.06 ±0.04 

0.03 ± 0.03 

0.11 ±0.02 

0.04 ± 0.04 

0.14 ±0.04 


Note. — Regions A and B contain flux from only the fraction of the source within the tangential caustic. 
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Fig. 2.— SED of the four lensed images and central AGN of the 
lens galaxy. The points are offset in wavelength for clarity. The 
shape of the AGN SED is distinct from that of the background 
source, indicating that the emission is from the lens galaxy rather 
than a central image. 


profiles, each defined as 


I(x,y) = A 


4 ur 


(i + q ) 2 


( 1 ) 


where (x, y ) are the coordinates along the semi-major 
and semi-minor axes from the center of the light distri¬ 
bution, A is the amplitude, q is the axis ratio, <f> (not 
shown) is the position angle of the light distribution, w 
is the size of the core, and s is the power-law slope. The 
total light profile is /total = h + h with structural pa¬ 
rameters (A 1 ,q 1 ,(j) 1 ,wi,si) for/t and (A 2 , q 2l fc, u> 2 , s 2 ) 
for I 2 . We use two power-law light profiles since a single 
one was inadequate with significant image residuals. The 
centroids of the two power-law light profiles are set to be 
the same, and the other profile parameters are indepen¬ 
dent. 

Table [3] lists the parameters for the two light profiles. 
There is one central component (/ 1 ) and a more diffuse 
component ( I 2 ). The core radius of the lens light distri¬ 
bution is therefore set by that of the central component, 
i.e., Rcore = W\. We present fits for both HST bands as a 
consistency check, although we primarily use the F160W 
results because it is the higher signal-to-noise image and 
should be a better tracer of the galaxy’s stellar mass 
distribution than the bluer F110W band. We also fit a 
double power-law plus point source (convolved with the 
point spread function) to account for possible low-level 
AGN emission. The double power-law plus point source 


WFC3/F110W WFC3/F160W 



- 2-1012 - 2-1012 
Aoc (arcsec) 


Fig. 3.— Top: HST/ WFC3 imaging of SDP.81 in bands F110W 
(left) and F160W (right) on a 0.065"/pixel scale. Middle: Best- 
fit models (double power-law plus point source) to the galaxy light 
distribution. Bottom: Residual image when our model is sub¬ 
tracted from the data. The cyan annulus indicates the region that 
was masked out during the light fitting. The best-fit models give 
an inner core radius of 0.15”, while the models without a point 
source give an inner core radius of 0.08 ,/ . All images are 5.655 ,/ on 
a side. 


fit has x 2 « 1500 in F110W and x 2 « 2300 in F160W, 
while the double power-law without a point source has 
X 2 ~ 1600 in F110W and x 2 ~ 2600 in F160W. 

We test a Sersic profile fit, but find that it is disfavored 
in the fitting region (x 2 ~ 1700 in F110W, x 2 ~ 2900 in 
F160W). We also test a double Sersic profile with linked 
centroids. The quality of the fit is nearly the same as the 
two power-law model, and the Sersic indi ces are consis¬ 
tent w ith those of the same profile fit by iTamura et al.l 
(|2015D to the F160W data. We prefer the two power- 
law model, as it is easier to determine a core radius ( w\) 
and more straightforward to associate to the mass model 
parameters for the composite model (Section 13. 3. 21) . 

We constrain w\ to be 0.08 ± 0.01" (« 356 pc) for 
the two power-law model, and 0.15 ± 0.01" (~ 668 pc) 
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TABLE 3 

HST Lens Light Fitting 


Parameter 

Posterior (Double Power-Law) 

Posterior (Double Power-Law + Point Source) 

F110W 

F160W 

FI 10W 

F160W 

Ai 

0 fiS'S - * -0 030 
u.ooo_o Q38 

0 452”*” 0 ' 011 
u.^oz _ 0 017 

U.OOZ_Q 078 

0 466 +0 ' 0081 
u.^oo_ 0 0082 

<71 

0.82 ± 0.01 

0.81 ± 0.01 

0.82 ± 0.01 

0.81 ± 0.01 

01 (°) a 

15 ±1 

12 ±1 

14 ±1 

12 ±1 

wi (") 

0.080 ± 0.003 

0.084 ± 0.001 

0.148 ±0.011 

0.148 ±0.004 

Si 

0.94 ± 0.02 

0.93 ±0.01 

1.04 ±0.05 

1.01 ± 0.01 

A 2 

3 O’*” 1 ' 1 
°* u -0.7 

0 4 + 1.4 

0.^_1 2 

6 7+4.3 
u ‘ 1 —2.4 

7 c+0.9 

0.9 

<72 

0.76 ±0.01 

0.78 ±0.01 

0.76 ±0.01 

0.78 ±0.01 

02 (°) a 

2 ± 2 

8 ± 1 

5 ± 1 

9 ± 1 

™2 (") 

1.25 ± 0.12 

1.54 ±0.12 

1.65 ±0.18 

1.93 ±0.05 

S2 

1.58 ± 0.14 

1.80 ±0.15 

1.70 ±0.20 

1.93 ± 0..05 


Note. — Uncertainties are statistical only. Reported values are medians, with errors corresponding to the 16th and 84th percentiles. 
a Angles measured East of North. 



Fig. 4.— High-resolution Band 7 continuum negative image. The 
beam size is the black ellipse in the bottom left. Open symbols are 
the positions used for the lens modeling for feature 1 (green trian¬ 
gles), 2 (magenta circles) and 3 (red diamonds). Filled symbols are 
the positions of the corresponding images predicted by our most 
probable power-law lens model for the R core = 0.08” case (the po¬ 
sitions from the R cor e = 0.15” model and composite model are very 
similar). The small filled symbols near the center are the modeled 
source positions for the same case. The critical curves (solid blue) 
and caustics (dashed orange) are overplotted. The image is 4.2” on 
a side. 

for the two power-law model plus a point source (where 
the uncertainties are obtained by adding in quadrature 
the statistical uncertainties and the systematic uncer¬ 
tainty from the difference in the two bands). We test 
both cases as priors on the core of the mass profile for 
the cored power-law model, as the compact central flux 
could be due to either AGN activity or stars. We also 
use the light distribution from these fits to constrain the 
stellar component of the co mposite model. We no te that 
an independent analysis by iTamura et al.l (120151) find a 
similar core radius (~ 0.15") from a fit to the same HST 
data. 


3.3. Lens Mass Modeling 


We try two models for the mass distribution of the lens 
galaxy: a cored elliptical power-law mass distribution, 
and a composite model with separate stellar and dark 
matter components. For each model, we sample the pos¬ 
terior probability distribution of the lens parameters r) 
given the observed data d. 

P{ri\d) oc L{d\rf)Q{ij). (2) 

The likelihood L is 


L(d\rj) = L pos (d|j 7 )L flux (d|rj), (3) 

where L pos is associated with the observed image posi¬ 
tions A, B, C and D, and Lfl ux is the non-detection of 
the central image E. Specifically, 


Lpos — ry 

4/, 


■ exp 


pos 


1 


N { 


e 

j =1 i=A,B,C,D 


| oobs _ ppred 


hJ 




(4) 

where Nf = 3is the number of features/clumps identi¬ 
fied in the background source, R°j S = (Xij s ,y°j s ) is the 
observed image position (listed in Table [T}, R?™ (rj) is 
the predicted/modeled image position, ay, is the uncer¬ 
tainty in the observed image position, and Z pos is the 
normalization given by 


N { 


-'pos 


= ( 27r ) 




n n 

j =1 i=A,B,C,D 


2 

a i,j 


(5) 


with N pos = 8 as the total number of image positions 
identified. The non-detection likelihood of the central 
image E in a CO map is 


fflux - 


V2^a C o 


exp 


/ r d (??) 2 

2<t co 


( 6 ) 


where aco is the background rms in the CO image and 
/r d (*7) is the predicted flux density of image E that 
is obtained by multiplying the observed integrated flux 
density of image D with the modeled magnification ra¬ 
tio iie/i^d- As discussed in Section I3TT1 we use image D 
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because it is the image of the entire source that is fur¬ 
thest from the critical curves, making it less susceptible 
to differential magnification effects. 

We calculate non-detection likelihoods for each of the 
CO J=5-4, CO J= 8-7, and CO J=10-9 maps using 
Equation then multiply them to get the final likeli¬ 
hood. In practice, we first sample L pos x Q via MCMC 
of chain length 5 x 10 5 , and then weight these samples 
by Eflux- We incorporate uncertainty in the flux density 
of image D by drawing from a Gaussian distribution set 
by the measured values and uncertainties (see Table [2]). 
We use the CO flux densities because there is no con¬ 
tamination from the AGN of the lens galaxy after the 
continuum subtraction. 

With our lens model, we map the observed positions of 
the three identified source features (Section Id.II) to the 
source plane. For each feature, we take the weighted av¬ 
erage of these mapped locations to be its position in the 
source plane. The mapped source positions are weighted 
by y/Jti/ai, where /q is the model magnification at the 
position of image i and eq is its associated uncertainty 
(Table [l]). To calculate the model magnification ratio 
He/Hd, we use these source positions and place circu¬ 
lar Gaussian profiles to create a mock extended source, 
consisting of three distinct clumps, that approximates 
the intrinsic source brightness distribution in the Band 
6 and 7 continuum th a t has bee n independ e ntly d eter- 
mined bv lRvbak ct all (12015H and lDve et al.l (12015t h We 
then calculate the magnification ratio he/hd by lens- 
ing this mock source with our model and calculating the 
relative flux of the D and E components. Additional 
tests show that our results are unaffected by moderate 
(~ 30 — 40%) uncertainties in the assumed source sizes. 

There is a satellite galaxy ~ 4.5" from the lens center, 
which is close enough that its influence may not be ade¬ 
quate ly described by external shear alone (|McCullv et al.l 
120141) . Its integrated light in F160W is ~ 1% that of the 
lens galaxy. We test models with this galaxy included 
and find that our results are unaffected. 

3.3.1. Cored Elliptical Power-Law Model 

We model the lens galaxy as a combination of a cored 
elliptical power-law mass distribution, a point mass rep¬ 
resenting the central SMBH, and external shear. The 
lens model parameters are given in Table [4] 

The prior Q{rf) in eq uation ([2]) on m ost of the pa¬ 
rameters are taken from iDve et all (120141 ) based on the 
HST/ WFC3 data (see Tabled]). For the prior on R corc , 
we test the two extreme cases determined from the lens 
light fitting (Section 13.21) . R CO re = 0.08 ± 0 . 01 " and 
Rcore = 0.15 ± 0.01". We assume that the SMBH is 
located at the center of the mass distribution and adopt 
a uniform prior on its logarithmic mass for both R CO re 
models. 

3.3.2. Composite Model 

We also fit a composite lens model that includes sepa¬ 
rate components for the stellar and dark matter contribu¬ 
tions, as well as a point mass for the SMBH, and external 
shear. We take the power-law fits to the stellar light dis¬ 
tribution from our best fit model to the HSTF160W data 
(Section 13.21) and use the parameters of those profiles to 
represent the stellar mass. We allow for Gaussian uncer¬ 
tainties on the profile parameters (Table [3]). Although 


the best-fit light model parameters are not exactly the 
same as the median values reported in Table [3] they are 
typically within ler. We use the stellar light parameters 
from the fit that included the central point source, as 
it was a better fit to the HST data. We represent the 
dark matter halo with a spherical NFW profile. From 
the results of iGavazzi et al.l (|2007l) . which are based on 
the S loan Lens ACS Survey (SLACS; iBolton et al.ll2006l 
I2008D sample of lenses with similar redshifts and velocity 
dispersions to SDP.81, we set a Gaussian prior on the 
scale radius of the NFW halo. In practice, this prior has 
little effect on the modeling since the scale radius is an 
order of magnitude larger than the Einstein radius of the 
lens. The priors on the position of the mass centroid and 
on the external shear are the same as those for the power- 
law models. Tabic |4] contains the model parameters that 
are not tied to the light distribution. 

The stellar mass-to-light ratio is assumed to be con¬ 
stant across the galaxy. We enforce a Gaussian prior on 
the total slope of the mass model at th e Einstein radius 
of T = —0.465 ± 0.03 (IDve et al.ll2014l) . T is defined to 
be ( 7 ' — l)/ 2 , where the three-dimensional mass density 
is p(r ) oc r -7 . 

4. RESULTS 

4.1. Lens Modeling Results 

We present the results of our lens mo del fitti ng in 
Table tD Our results are consistent with iRvbak et all 
(I2015D . who fit directly to the data in visibility space. 
Comparing the “Prior” and “Posterior” parameter val¬ 
ues, most parameters are better constrained by the multi¬ 
ple image positions identified from the ALMA data. Our 
models predict the central image to be near (< 0.05") the 
mass centroid (see Table 0]). Since it is not detected, the 
sampling is relatively insensitive to R CO re for the power- 
law models, so the posterior of R CO re is essentially set by 
the prior determined from the HST/ WFC3 data. 

Most parameters for the two different R CO re priors for 
the cored power-law model are very similar since the dif¬ 
ferent core radii are not affecting the overall fit to images 
A-D. In Figure^ we show the predicted image positions 
of the four features in A, B, C and D of the most proba¬ 
ble model with the R CO re = 0.08" prior. For both models, 
we can reproduce the observed image positions with an 
rms offset between the observed and predicted image po¬ 
sitions of ~ 0.015". 

The composite model provides an acceptable fit (rms 
image offset ~ 0.031"), but it is not as good as the power- 
law models. The predicted stellar mass-to-light (M/L) 
ratio is l.l^o i M 0 Lq 1 at a rest-frame wavelength of 
~ 1.2 pm, which is consiste nt with the M/L at rest- 
frame ~ 1 pm determined by iTamura et all (120151 ) that 
is needed to normalize the stellar light and mass profiles 
in their lens model. The dark matter fraction within the 
Einstein radius is ~ 50%, which is consistent with and on 
the high side of various lensing and/or dynamical studies 
of similar lens g alaxies at a similar fraction of the effec¬ 
tive radius (e.g., Barnabe et al.ll 20 TTl : iQguri et all 12014 
iSonnenfeld et aljl2015l). 

We show the convergence profiles for the three mod¬ 
els in Figure [5] Since the central image magnification is 
inversely proportional to the surface density in the inner¬ 
most regions, we expect, from these results, to have the 
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TABLE 4 

Lens Model Parameters 



Cored Power Law Model 


Parameter 

Prior 3, 

Posterior (0.08" R C ore prior) Posterior (0.15 // R C ore prior) 

a (J2000) 

Gaussian; 09:03:11.572 (= 0.0) b ±0.005" 

09:03:11.572 (-0.002) b + 0 001 " 09:03:11.572 (-0.002) b + 0 001 " 

-0.005" v ' -0.005" 

5 (J2000) 

Gaussian; +00:39:06.54 (= 0.0) b ±0.005" 

+00:39:06.53 ( 0.002) b + 0 005 '' +00:39:06.53 ( 0.002) b + 0 005 '' 

v ' -0.005" v ' -0.005" 

0 E (") 

Gaussian; 1.56 =b 0.12 

1.60 ± 0.01 

1.60 ± 0.01 

r c 

Gaussian; 0.465 ± 0.03 

0.47 ±0.03 

0.47 ±0.03 

b/a 

Gaussian; 0.79 =b 0.04 

0.82 ± 0.03 

0.82 ± 0.03 

e (°) d 

Gaussian; 10.8 zb 6.9 

13 ±4 

13 ±3 

Rcore (") 

Gaussian; (0.08; 0.15) ± 0.01 

0.08 ± 0.01 

0.15 ±0.01 

Rcore (P C ) G 


356 ± 45 

668 ± 45 

'Text 

Gaussian; 0.04 =b 0.02 

0.03 ±0.01 

0.02 ± 0.01 

0-y (°) d ’ f 

Uniform; [— 00 , 00 ] 

00 

1 + 

-4 O 

74+ * 1 2 
'^-21 

log(M B H)(MQ) 

Uniform; [6.5,9.5] 

(see Figure RTb 

(see Figure l6l) 


Composite Model 


Parameter 

Prior 

Posterior 


Stellar M/L (Mq L+s 

Uniform; [0,oo] 

i-i+l 


NFW k s 

Uniform; [0, oo] 

0-078++ 


NFW r s (") 

Gaussian; 18.6 =b 2.6 

i8.3ll;5 


'Text 

Gaussian; 0.04 d= 0.02 

0.06 ± 0.01 


e-y (°) d ’ f 

Uniform; [—oo, oo] 

-set) 


log(M B H)(MQ) 

Uniform; [6.5,9.5] 

(see Figure [ 6 } 



Parameter 


Predicted Central Image Position 
Posterior (0.08 ,/ Rcore prior) Posterior (0.15" R CI 


prior) 


Posterior (Composite) 


77 “ 


a (J2000) 
<5 (J2000) 


09:03:11.573 (-0.020) b 


+ +0.007" 
-0.007" 

+00:39:06.54 ( 0.010) b + 0 006 '' 
V ’ -0.005" 


09:03:11.574 (-0.037) b 
+00:39:06.55 ( 0.016) b 


+0.008' 

-0.009' 

+0.007' 

-0.006 


n 


09:03:11.573 (-0.014) b + 0 006 
v ' -0.006 

+00:39:06.54 ( 0.010) b + 0 005 '' 
K ' -0.005" 


Note. — Reported values are medians, with errors corresponding to the 16th and 84th percentiles. 
a Model priors from lDve et all (2014 1 based on HST/ WFC3 data. Core radius prior is from our fits to the WFC3 lens light distribution. 
For Gaussian priors, the center and width are given. For uniform priors, the range is given. 
b Values in parentheses indicate offset in arcseconds relative to the mass centroid prior. 
c r = ( 7 ' — l)/ 2 , where the three-dimensional mass density is p(r) = r —7 . 
d Angles measured East of North. 
e Physical Rcore scale under our assumed cosmology. 

i 6 1 = 0° corresponds to shearing along North-South direction (i.e. external mass distributions East or West from the lens system). 
s Stellar mass-to-light ratio is given for the WFC3/F160W band in solar units at roughly rest-frame 1.2 //rri (corresponding to observed 
F160W at z = 0.2999). There is no direct prior on the M/L, but the value is set by enforcing a Gaussian prior on the slope of the combined 
mass profile of T = —0.465 ± 0.03 IIDve et al.1120141 ). 


strongest constraints on the SMBH mass from the cored 
power-law model with R CO re = 0.15" and the weakest 
constraints from the composite model. 

4.2. Limits on SMBH Mass 

The non-detection likelihood in equation [6] places lim¬ 
its on the mass of the SMBH. Figure [G] shows the relative 
posterior probability density of the SMBH mass. For the 
power-law model with R cor e = 0.08" and the composite 
model, the lens galaxy has a sufficiently high central den¬ 
sity that the predicted central image is substantially de- 
magnified below the flux limit, so the non-detection is un¬ 
able to add useful information to the SMBH mass. How¬ 
ever, for the power-law model with R CO re = 0.15", the 
probability density rises sharply at Ic^Mbh/Mq) ps 8.5. 
This is the mass at which the central image demagnifica¬ 
tion transitions from being controlled mainly by the mass 
density of the cored power-law mass distribution to be¬ 
ing controlled mainly by the SMBH. The hatched region 


in Figure [G] shows the 1 -a range of SM BH masses pre¬ 
dicted by the MsH-Mbuige relation of iKormendv fe Ha 
(12013T) for SDP.81, assu ming the stellar mass determined 
bv lNegrello et al.1 ( 2014 ) scaled for a Salpeter initial mass 
function. For R cor e = 0.15", the non-detection of the cen¬ 
tral image can tighten the range of likely SMBH masses. 
The solid grey region sho ws the same relation assuming 
MeH/Mbuige ex (1+z) 1 ’ 96 dBennert et al.|[20TH) . although 
the redshift evolution of the MeH-Mbuige relation is fairly 
uncertain from current observations. 

4.3. Constraints from Deeper Observations 

With deeper observations, tighter constraints can be 
made if the central image is still undetected. In the lower 
panels of Figured we show results for a hypothetical ob¬ 
servation of SDP.81 with a S/N ratio three times higher 
than the current CO data. The ALMA observations 
are starting to be able to constrain the R CO re = 0.08" 
power-law model and the composite model, and the con- 
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Fig. 5. — Convergence profiles for the cored power-law model 
with Rcore = 0.08 // (blue solid line), Rcore = 0.15 // (red dashed 
line), and composite (grey dash-dotted line) models. We expect 
models with a lower central surface density to provide stronger 
constraints on the SMBH mass. 

straints for the R CO re = 0.15 r/ power-law model much 
more strongly rule out log(MBH/M Q ) < 8.5. The tran¬ 
sition SMBH mass for the smaller core radius and the 
composite model is lower, log(M B H/M Q ) ~ 8.0 — 8.2, 
indicating that a smaller black hole mass is needed to 
demagnify the central image. 

In Table SJ we also show the predicted location of the 
central image. These positions and their uncertainties 
are calculated by taking the weighted median and 16 — 
84% quantiles, but only for models in which the central 
image is detected. The predicted position of a central 
image in SDP.81 is very close to the lens center (< 0.05") 
for all models. If it were detected, it may be possible to 
use the position of the central image as an additional 
constraint on the lens model. Such a constraint would 
be unlikely to impact the large-scale model parameters, 
but could place additional constraints on R CO re and the 
SMBH mass. In principle, it could also help to constrain 
the lens centroid or a possible offset between the mass 
centroid and the SMBH position, which we have assumed 
to be coincident in our model. Since it is undetected 
in the case of SDP.81, a further investigation of these 
possible constraints from the central image position is 
beyond the scope of this paper. 

Nonetheless, we explore specifically the constraints on 
the SMBH mass for the models we have considered, us¬ 
ing the flux measurement of a hypothetical detection of 
the central image. Such a detection could place a direct 
constraint on the SMBH mass rather than just a lower 
limit. In Figure [71 we show the SMBH mass constraints 
for our three models assuming a hypothetical detection 
of a central image with a flux of 0.3 and 0.1 times the 
background rms level of the current data at a S/N of 
3. The position of the central image is not used as a 
constraint. Depending on the assumptions of the lens 
model, a detection of the central image can either di¬ 
rectly constrain the SMBH mass to within a range (e.g., 
the power-law R cor e = 0.15" model), or provide an upper 
limit (e.g., the power-law R cor e = 0.08" model and com¬ 
posite model when the central image flux is 0.3 times the 


background rms level). A detection sets an upper limit 
because larger SMBH masses would demagnify the im¬ 
age below the observed flux, while the lower limit (in the 
cases of a direct constraint) comes from demagnification 
due to the lens galaxy profile itself. This highlights the 
importance of complementary observations to constrain 
the central mass distribution of the lens galaxy, such as 
high-resolution imaging to constrain the stellar surface 
mass density, or dynamics. 

5. CONCLUSIONS 

Combining high-resolution long-baseline ALMA data 
and archival HST imaging, we model the gravitational 
lens SDP.81 and investigate the prospects for detecting 
the demagnified central image. We detect compact emis¬ 
sion region in the central region of the lens, but its SED 
indicates that it arises from low-level AGN activity in the 
lens galaxy. There is no evidence for a central image in 
any molecular lines. Using the positions of three distinct 
features in the source, we model the mass distribution of 
the lens galaxy with both a cored power-law model and 
a composite model including separate stellar and dark 
matter components. 

Based on the non-detection of the central image in the 
CO maps, we are unable to constrain the SMBH mass 
for the Rcore = 0.08" power-law model and the compos¬ 
ite model, but the R CO re = 0.15" model shows a prefer¬ 
ence for a SMBH mass of log(MBH/M 0 ) > 8.5. Deeper 
ALMA observations can strengthen this constraint and 
place limits on the SMBH mass for models with higher 
central density. A deeper observation that is able to de¬ 
tect the central image can place a direct constraint on 
the SMBH mass, or an upper limit, depending on the 
assumed lens galaxy profile. These hypothetical con¬ 
straints highlight the importance of complementary ob¬ 
servations to constrain the central surface mass density 
of lens galaxies for this purpose. 

Future ALMA observations of strongly-lensed galaxies, 
particularly two-image lenses whose central images are 
less-strongly demagnified, may yield better opportunities 
for detection. Shorter baseline observations will also be 
more sensitive and could improve the chance of observing 
a central image in the molecular lines. Such a detection 
would place interesting constraints on the central mass 
distributions of galaxies at cosmological distances and 
serve as an independent probe of their SMBH mass to 
better understand the galaxy-black hole connection and 
its evolution over cosmic time. 
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PL, R core = 0.08 PL, R core = 0.15 Composite 



Fig. 6. — Constraints on the SMBH mass from the non-detection of the central image in CO. We show constraints for the cored power-law 
models with Rcore = 0.08 // (upper left) prior, R C ore = 0.15" (upper middle) prior, and the composite model (upper right). The red line 
shows the prior probability density, which is uniform in log(MBH/M0). The blue line shows the relative posterior probability density 
when the non-detection is taken into acco unt. Th e normalization is arbitrary. The hatched band sh ows the l-q range of S MBH masses for 
SDP.81 from the MBH-Mbuige relation of IKormendv fc Hoi (1201311 , assuming the stellar mass from [Negrello et al.l (120141 1, while the solid 
grey band is the same quantity accounting for redshift evolution of MgH/Mbuige oc (1 + 2) 1,96 JBennert et al.|[2011l) . The non-detection 
of the central image indicates a preference for log(MBH/M©) > 8.5 for the Rcore = 0.15 // case, which can improve the constraint from 
the MBH - Mbuige relation alone. The R C ore = 0.08" and composite models do not provide additional constraints on the SMBH mass, as 
the central density is too high. The lower panels show the results for an observation with a S/N ratio three times higher in the CO lines. 
With these hypothetical constraints, we can strongly rule out log(MBH/M0) < 8.5 for R C ore = 0.15". We start to see a preference for 
log(MBH/M0) > 8.2 when R C ore = 0.08" and log(MBH/M0) > 8.0 for the composite model. 

NRC (Canada) and NSC and ASIAA (Taiwan), in coop- eration with the Republic of Chile. The Joint ALMA Ob¬ 
servatory is operated by ESO, AUI/NRAO and NAOJ. 
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Fig. 7. — Constraints on the SMBH mass in SDP.81 from a hypothetical detection of the central image for the cored power-law models 
with Rcore = 0.08 // (left) prior, R C ore = 0.15 // (middle) prior, and the composite model (right). We show constraints from a detection of a 
central image with a flux of 0.3 times the background rms level of the current data (blue) and 0.1 times the background rms level of the 
current data (red), both with a S/N of 3. The normalization is arbitrary. Depending on the assumptions of the lens model, a detection of 
the central image can either directly constrain the SMBH mass to within a range (e.g., the power-law R C ore = 0.15 ,/ case), or provide an 
upper limit (e.g., the power-law Rcore = 0.08" case when the central image flux is 0.3 times the background rms level). 
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